# 离差标准化处理
min.max.norm <- function(x){
  ((x-min(x))/(max(x)-min(x)))
}

# 读数据到x
x=scale(data)
R=round(cor(x),3)
eigen(R) # 取3个因子
a=factanal(x,factors = 2, rotation = "varimax") # 极大似然法做因子分析
b=factanal(x,factors = 2, scores = "regression") # 回归法计算因子得分
print(b$scores)